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Abstract 

04 , The complete model of production of baryons in expanding primordial 

■ symmetrical hot Universe is constructed in terms of general relativistic 
D ' kinetic theory. 

Q' 

1 Introduction 

An important example of local thermodynamical equilibrium (LTE) breaking 
by massive particles during the cosmological expanding is process of generation 
^ ' of baryons (baryogenesis) in primordial baryon - symmetrical Universe. On pQ, 

Vh^, [5] there was promoted a hypothesis that observing baryon skewness of Universe 

5b = «£ „ io-9 (1) 

> ' 

■ is caused by CP- non - invariant processes, which break the conservation of 
y—{ , baryon's charge. An example of such process is decay of superheavy X- bosons, 
\Q ■ to which Grand Unified Theories (GUT's) adduce^ 

in 

(Nj i X^q + q- X^q + l, (2) 

where q - quark's symbol, I - lepton's symbol. In 3 , however is shown that if 
given processes are run in LTE's conditions then even the existence of - non- 
invariance and the breaking of baryon charge's conservation is insufficiently for 
appearance of baryon's excess over anti-baryons. For production of the baryon's 
charge in primordial baryon-symmetrical Universe besides stated facts require 
the breaking of LTE in reactions of such type ^ and withdrawal of baryons 
from statical equilibrium's conditions. Such possibility: 

Tx > t (3) 

(t x - time of half-decay of X - boson, t - cosmological time ) is realized under 
condition 0]: 

m x > a x m p iVN, (4) 



X 



(mx - mass of X-boson, m p i = y § - Planck's mass, N - number of degrees 
of freedom). Condition Q stringently limits from below allowable values of X- 
bosons masses. 



1 Here and further a means antiparticle a. 
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In Ref. [3], [5] (see too [7], [5] and etc.) there were made estimates of baryon 
skewness of Universe with the assumption of LTE's breaking. For attitude of 
baryons number's density to density of total entropy, 5, there was obtained: 




n b 45C(3) iVx 
S ~ 4tt 4 N 



Ar, 



(5) 



where ((x) - £ - Rieman's function, Nx - the number of types of superheavy X- 
bosons, Ar - difference of comparative probabilities of decays in channels 



which appears in consequence of breaking of CP - invariance. Further in papers 
[5]-[TT] on basis of relativistic kinetic theory there were made numerical accounts 
of attitude ns/S, which basically confirmed foregoing estimates. In cited earlier 
papers on the basis of made calculations there was found the inferior limit of X 
- boson's mass: 



which created considerable difficulties for standard SU (5) theory. 

Previously mentioned papers however have considerable deficiency - instead 
of direct solution of kinetic equations for X - bosons under conditions of con- 
siderable LTE's breaking boson functions of distribution are described by quasi- 
cquilibrium (Bolzman's) distributions, which parameters are defined from hy- 
drodynamical equations for moments of distribution's function, i.e., in this pa- 
pers actually is used hydrodynamical method of Grad. According to results of 
relativistic kinetic theory (see for example, [H]-[TI]); global thermodynamical 
equilibrium in homogeneous isotropic expanding plasma is reached only in ex- 
treme non-relativistic limit, or in extreme ultrarelativistic limit. In range of 
intermediate energies of particles and on conditions that LTE is breaking, the 
distribution of particles is not approximated by equilibrium distribution. In 
Ref. [UJ-[IB] of one of authors in the framework of kinetic theory there was 
found non-equilibrium function of distribution of X - bosons and shown that 
this function may essentially differs from equilibrium function. Therefore results 
obtained in [9j-|lljare correct if strong inequality Q is applied, but needed a 
correction in other areas. The estimate ([3]) although seems to appear sufficiently 
accurate, but doesn't describe such situations, when condition is not realized 
while width of covering of experimental and theoretical implications £5 permits 
possibility when baryon charge is produce in conditions, that are less favorable 
than Q. 

The answer to the question, what kind will be magnitude 63 at mx < 
oexmpiVN, could be given only by detailed kinetic analysis. On the assumption 
of confidence to the hypothesis of primordial baryon- and charge- symmetrical 
Universe, such analysis from the other side allows more definitely indicate range 
of possible implications of fundamental constants of Grand Unified Theories. In 
the middle of 80-s - beginning of 90-s by one of authors was formulated the 
kinetic model of describing of symmetry's breaking processes and there were 



X -t q + l X ->• q + l 



m x > 10 lb Gev, 



(6) 
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obtained some estimations, following from this model (Ref. [T7J-Q35])- Specif- 
ically on basis of obtained estimations there was formulated a hypothesis that 
accounting of kinetic of baryogenesis process will permit to reduce the low esti- 
mation of masses of superheavy bosons to quantity of order 5 • 10 14 Gev. However 
there wasn't executed detailed analysis of this model in this papers and on ac- 
count of external, researches of this problem weren't complete. The purpose of 
this paper is exactly carrying out of more detailed analysis of the kinetic model 
of baryogenesis and construction of the numerous model of considered events. 
Let us notice that though this paper contains results for concrete model of in- 
teractions based on the minimal SU(5) symmetry, the generalization of these 
results for other analogical field theories isn't hard - it is reduce to arithmetical 
recalculation of corresponding coefficients. 

2 The algebra of interactions 

Let us consider for example standard SU(b) model of interactions (see for ex- 
ample 8 ). There are involve 12 vectorial calibrating bosons which represent 
themselves 2 charged color triplets in this model: 

r -yi y-i y-i y-i \ 

t^V-4/3' ^,4/3' ■ A ti,-l/3>- A «,l/3.r 

- i - color index (red, green, blue), the subscript - electrical charge - and 12 
scalar Higgs's bosons: 

r y-i y-i y-i y-i t 

t^s,-l/3' ^s,l/3' ^8,-4/3) ^-s,4/3J- 

In more common models the number of vectorial bosons reduces to 24, - and 
2 more triplets add at the same time 

f yi yi . yi yi \ 

X^v',-1/31 ^-ti',1/3' ^u',2/3' ^v',-2/3J i 

and the number of Higgs's bosons reduces to 30, and at that add 3 SU(3) - 
triplets, which are included in 3 SU (2) triplets: 

f yi y-i . -yi y-i . yi yi \ 

VV,-l/3> ^s',1/3' ^-s',2/3' ^s'-2/3' ^s',-4/3' ^s',4/3/- 

Superheavy X-bosons are often calling X, Y, Z-bosons according to their charges: 
—4/3, —1/3, 2/3. The electric charge's conservation laws allow to run only these 
reactions of decay /creation of these bosons: 

^-4/3 Qa + <la\ -X"-4/3 q K + ^e! 

X-l/3 ^ Qa + 9k5 X-l/3 ^ Ik + lv) -X"-l/3 ^ Qa + l e ] (7) 
-^2/3 ^ Qk + 9k 5 X 2 /3 Qa + lis 

where a - index of apoquark (a = u, c, t - quarks with charge 2/3), k - katoquark 
{k = d, s,b quarks with charge —1/3), e - symbol of charged lepton (e = e, /i, r_, 



3 



charge=-l), v - symbol of neutral lepton (y = v e , u^, v T ). In reactions ([7} color 
charge is also conserving so therefore in reactions of decay of X-boson and in 
reaction of quark annihilation antiquarks of different colors, additional to X 
-boson's color, take part in. For example: 

*-4/3 = d r + e- = «f + uf ; 

where R,G, B - indications of red, green, blue colors; corresponding overlined 
indexes R,G, B - anticolors, which are equal to sums of corresponding addi- 
tional colors, l 7 r - indications of left and right particles. We may suppose, as 
it often doing that in right parts of reactions (J7]) particles only with the same 
charm take part in. In that case in standard variant SU(5) there will be 90 cou- 
ples of reactions of type (0, and in expended variant of theory - 198 couplefl 
Then, the number of degrees of freedom in standard variant of theory is equal: 

N = J2(2s+l) + l + = 185 

B F 

- 24 vectorial bosons (s — 1), - 29 scalar bosons (s = 0), 36 quarks (s = 1/2), 
12 leptons (s = 1/2). In expanded variants SU(5) this number is even greater. 

3 Kinetic equations 

Let us for convenience express averaged on spin states invariant elements of tran- 
sition matrix of decays of superheavy bosons with the help of nondimensional 
numbers (r^, fi): 

\M x ^ q if = \s 2 r e ; \M x ^ q if = \s 2 r v ; 

\ M x-Hfif = 1^5 \ M x^gL T = §* 2 *V5 ( 8 ) 

\M x ^ m \ 2 = \s 2 (\ - r); \M x ^ qq \ 2 = ±s 2 (l - f), 

(r = r e + r„, f = f e + f„), which should be also provide by indexes of electrical 
and color charges as well as by charms of quarks and leptons. 

Through CPT - invariance averaged matric elements of reverse transition 
(annihilation) are equal: 

\M q i^ x \ 2 = \s 2 f e - \M ql ^ x \ 2 = \s 2 f u - 

\M q j^ x \ 2 = \s 2 r e ; \M ql ^ x \ 2 = \s 2 r„; (9) 

\M qq ^ x \ 2 = \s 2 (l - f); \M qq ^x\ 2 = h 2 {l - r). 

2 With account of decay and creating of superheavy antibosons. 
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Invariant matric elements of double-particle decay's probability are constant. 
Calculation of total probability of decay in one-loop approximation in the con- 
text of, for example, standard SU (5) model gives: 

W x ^ ql ^^(2irfm 2 x 5^{P F ~Pj), 

where P-total 4-momentums of initial and final conditions a ~ 1/45. In that 
case: 

A,C Z ' 

where summation is carries over all charms and colors of particles participating 
in decay of given type of X -boson. 

In further however we shall not define concretely s 2 , supposing only: 

s 2 = Const. (11) 

Besides first order reactions ([7]) reactions of higher order, running with partici- 
pation of superheavy bosons and breaking CP-invariance and also conservation 
laws of baryon and lepton charges also are possible. For example 

qq' 4 ql; ql 4 qq' . (12) 

Let us consider further metric of homogenous isotropic Universe d: 

ds 2 = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ). (13) 

Let further p 1 - 4-momentum of a particle, corresponding to relation of nor- 
malization: 

(P,P) = (P 4 ) 2 - P 2 = m2 E — \Jm 2 +p 2 , (14) 

where p 2 = a 2 (t)((p 1 ) 2 + (p 2 ) 2 + (p 3 ) 2 ) - square of kinematic momentum of 
particle, E = p 4 - kinetic energy of particles. 

In Ref. |15) was shown that relativistic kinetic equations are asymptotically 
conformally invariant in ultrarelativistic limit: 

<P 2 > /, c x 

— — > oo (15) 

m 

to within 0(m 2 /p 2 ). Thus in ultrarelativistic limit (fl"5"j) when energies of all par- 
ticles participating in reactions lot more than their rest masses, kinetic equations 
in Friedman's metric after corresponding conformal conversion coincide with ki- 
netic equations in Minkovsky space, therefore homogenous isotropic expanding 
of ultrarelativistic plasma doesn't lead out from condition of global thermody- 
namical equilibrium. Further in opposite non-relativistic limit: 

^ - 0, (16) 
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Here and further we choose system of units h = c = G = 1 . 
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when kinetic energies of all particles participating in process far less than their 
rest masses, global thermodynamical equilibrium is again restores in plasma. 

mm- 

Let's introduce according to [TS] conformal momentums of particles P, Q, 
which are integrals of motion relative to metric (| 13[) according to formula: 

P = a(t)p (= V(Pi) 2 + (P2) 2 + (p 3 ) 2 ), (17) 
conformal kinetic energy of particles E = a(i)E: 

E = yja 2 (t)m 2 +P 2 (18) 

and corresponding macroscopic conformal parameters of plasma , which con- 
serve constant at the ultrarelativistic stage of expanding: 

T = a(t)T — temperature; (19) 

Af = a 3 (t)n— number density of particles; (20) 

£ = a A {t)e — energy density. (21) 

Kinetic equations for plasma in expanding homogenous Universe with par- 
ticipation of X -bosons we will write in symbolic form [19] : 



+ (other interactions), (23) 

where f a (P,i) - phase-space densities, F - index of fermions, P, Q - conformal 
momentums of X -bosons and fermions, correspondingly. Expression "other 
interactions" in (|23l) means integral of collisions for all other interactions which 
don't violate conservation of baryon and lepton charges: quark-gluon, lepton 
and quark decays of W- X-bosons, annihilation etc. All these reactions run in 
conditions of LTE and lead to establishment of thermal equilibrium of quarks 
and leptons. Integrals of interactions of particles we will write in relativistic- 
invariant form |20) : 

I a (x,P a ) = 

- E / 5(4) ( P / - P i)( Z n W if ~ ZifW if )tfdK, (24) 

f,i 

where summation is carrying out by all initial, i, and final, /, conditions of 
particles, integration (dsn 1 ) is carrying out by all four-dimensional momentum 
spaces, except space of particle of a- sort, W/i Wif - corresponding invariant 
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scattering matrixes which relate with invariant amplitudes of scattering Mi / by 
formulas |21) : 



W fi = (2Tr) 4 \M fl \ 2 ■ 2-( Ni+N ^- (25) 

- Ni Nf - numbers of particles in initial and final conditions. In formulas (|40p 
introduced products of initial and final conditions: 

^/i=n//iK i± /«)i 

Z if =Hf i l[{l±ff), (26) 

i f 

where symbol "— " relates to fermions, "+" - to bosons. Invariant elements of 
momentum space volumes are equal: 

d7la= (2^)3 *« P » **•)-—)> 

where p a - statistical degeneracy factor. 

Further, since reactions ([7]) . (|12[) run under very high temperatures T > 
10 14 Gev, it is possible with very high accuracy grade to mean that all fermions 
are ultrarelativistic, therefore: 

f Fa = [exp(-A Q + Q/T) + If 1 , (27) 

where conformal temperature T — aT, is constant in ultrarelativistic equilib- 
rium plasma, and adduced chemical potentials 

> _ 

are satisfy to series of conditions of chemical equilibrium, corresponding to reac- 
tions executing in conditions of LTE. Combining algebraic equations correspond- 
ing to these conditions we may come to such conclusion: chemical potentials of 
each type of quark aren't depend on their own color and charm. In that way 
it remains only 4 independent chemical potentials which for simplicity sake, 
following [8], we will denote by symbols of particles themselves - u, d, e, v: 

u = \ qa ] d = X qK ; e = A; e ; v = \ iv , (28) 

for them at that as for ultrarelativistic particles in conditions of LTE applies 
known condition of antisymmetry: 

A a = -X a . (29) 

In situation when rest mass of neutral leptons is equal to zero, their chemical 
potential is also equal to zero: v = 0, as chemical potential of massless parti- 
cles being in thermodynamical equilibrium Q In that case remain only three 
independent variables: u,d,e. 

4 see for example, 22 . 
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3.1 Kinematic equations for fermions 



At completion of chemical equilibrium conditions the last component in right 
part (p?3"|) becomes a zero. The second component in right part (j2"3"|) describes 
all reactions executing with violence of CP-invariance, in which quarks, leptons 
and virtual X-bosons participate. Since these integrals completely defined by 
equilibrium functions of distribution of quarks, leptons and other particles, then 
using functional Bolzman's equations |23] 

z fi -Zif = o => E A_ E A ' = ' ( 3 °) 

/ 

Fermion integrals of interaction can be write down in the form of: 

E /fF'-k.. = ~E J dTT F dn F dTr F ,f F f F , x 

x J\[d-K k {l± f k ){W FF ,^...-W FF ,^J. (31) 

There is no need to concretize probabilities of many-particle conversions W 
W for calculation of these integrals - it is enough to use optical theorem which 
is consequence of unitarity of S- matrix, [5] , [21] : 

e / n^ 1 ± - = °. ( 32 ) 

k 

where summation is carrying out by all final conditions of reactions FF' 
So in consequence of optical theorem we will obtain equality: 



^ FF - W x ^ FF ,)dir k + 

X,F> 

E'(l ± fk)(W FF ^... - W FF ,^Jd7r k = 0. (33) 

k 

(34) 

Using this result in equations (I23I) and integrating them by momentum space 
of fermions we will receive equations which in further will be used for definition 
of chemical potentials of fermions: 

dJV F . . - 



= a(t) J Wx->FF'dnxdTr F dTr F i x 



dt 

X,F' 



[f x (l - f F )(l - f F ,) - (1 + fx)f F f F '] ■ (35) 

Last equations can be still simplified by accounting of correlations which are 
correct to equilibrium fermion distribution functions (|27[) : 



J dlT F dTT F 'S^(px -Pf -PF>)fFfF> 
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2 3 7T 5 

where incorporated notations | 



= ^Efx(F + F')P(F,F'); (36) 
d-K F dTT F '5^\px -Pf -Pf')(1 - - /f') 

' [l + f x (F + F')]f3(F,F'), (37) 



f x (F + F') = [exp (-X F - Xpi + E x /T) - If 1 
P(F,F>) = p(F',F) = p(F) + P(F'); 

o (F] = T , i + cxp(-A F+P+ /r) i 

Pl j p l + exp(-A f +p_/T) 2' 



(38) 



p ± = -(£±p); S = £ x = ^m2 f +p 2 . (39) 

Function (3(F, i 7 ") is statistical factor which in Bolzman's statistics (A — > oo) is 
equal to 1. Arguments P and t of functions f x j3 are dropped for short. In 
necessary situations we will write explicitly: 

f° x (x;F,t); 0(x,F,t). 

Subject to (|36 ]) . ([37| equations (|35|) can be essentially simplified: 

M X ^ FF ,\ 2 x 



dJ ^F u\X^' PX 



= «(*)£': 



dt >j£ (27T) 

(40) 

Summation in ()40|) is carrying out by all X-bosons, moreover in quark-quark 
channels sum of colors of two quarks is complementary to color of X-boson, 
thus to each type of X-boson are correspond 2 addends in sum. 

3.2 Kinetic equations for X-bosons 

For X-bosons optical theorem assumes the form: 

J2 dir F dn F ,(l ~ f F )(l - f F >)(W x ^ FF , - W X ^ FF ,) = (41) 

F,F' 

and thus boson equations can be wrote in the form: 

x+r dt - 



in 



-^Y,\ M x^PF'\ 2 mF')[fx-f x (F + F') 

77 F,F' 



(42) 
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Let's draw attention to that important and strict fact, that kinematic equations 
for X-bosons become linear differential equations. Equations for antiparticlcs 
receive from (1401) (1421) by inversion of overlined indexes of particles - antiparti- 
cles. 

If assume that spontaneous violence of SU (5)-symmetry happened at very 
early stages of cosmological evolution when X-bosons were still ultrarelativistic, 
then mentioned system of kinetic equations (|40"]) . (|42|) is necessary to solve with 
initial conditions corresponding to initial global thermodynamic equilibrium and 
initial baryon and lepton symmetry: 

Aa| t=o = 0; /x(P,0) = / o o (0) = [exp(P/T o )-ir\ (43) 

where: 

T ° = T(0 >=Usv) ■ (44 > 

From (|4"2"| it is obvious that boson functions can be found in quadratures 
[2], [TB]. These equations in standard mathematical notation look like: 

fi + f i A i (t) = Y i (t), (45) 

point here and further means time derivative, i, k - now indexes of X-bosons, 
and vectors A Y define by correspondences: 

4(P,t) = —-= 4§=_^|M Xi ^| 2 /3(F,F'); (46) 



Yi (P,t) = ° W x 

J2 \M Xi ^ F<F ,\ 2 p(F, F')f x (F + F'), (47) 



- summation in (|46|) , (|47|) is carrying out by all channels of reactions over which 
given X-boson decays. 

Solving equations flU} with initial conditions (j43| . (j44|) . we will receive so- 
lutions for boson functions in quadratures: 



fx(P,t) = /x(P,0)exp(- / MP,t')dt') + 
K o ' 

cxp(- / Ai(P, t')dt'\ J ex P (/ Ai(F, t")dt")Yi(P, t")dt". 



Substitution of obtained solution for boson distribution functions in ki- 
netic equations for fermions (|40| leads to closed system of nonlinear integro- 
differential equations relative to chemical potentials. 
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3.3 Conservation laws 

In complete kinetic model of Universe it is necessary to add to equations (HQ)) , 
(j4"2")l another equations which define evolution of temperature and scale factor. 
First of these equations receives by integration of kinetic equations (j2"2")l . (j2"3"j) 
with weight Ei and following summation |23j : 

£ = a^mx(/Cx +IC X ), (49) 
x 

where £ - summary conformal energy density: 

7r 2 ivr 4 ^ p x 



30 ^ 2tt 2 



ES fxr 2 Ja^m 2 x +^dV, (50) 



2^ 2 7 ^c?m x + P 2 1 V ; 



where iV' - number of degrees of freedom of 5J7(5)-model without accounting 
of X-bosons: 

N' = N-N X , 

i.e., in standard SU(5)-modcl N' = 185 — 53 = 132. And equation that defines 
evolution of scale factor is Einstein's equation in which we can put k — 0, taking 
into account that we consider very early Universe t — > 0: 

a 2 = f£. (52) 

Equations (gDJ), (02J, (@SJ>, (EH jointly with definitions ((50]), and ini- 
tial conditions (j43|) . (J44J) are the complete system of equations which describe 
dynamics of baryogenesis. Let us consider some strict implications of these 
equations. 



1. Let us suppose that X-bosons lie in thermal equilibrium, i.e., ignoring kinetic 
equations for X-bosons (|4*2"j) . we will put in (14*U1) 



fx = f x {\x) 

where X x = —A — A'. Then exact solution of equations (I4U1) . satisfying to initial 
conditions (|4"3")l will be: 

N a = Const; X a = 0, (53) 

- i.e., in conditions of statistical equilibrium of X-bosons baryogenesis doesn't 
execute. 



exp 



-A x + 



U a 2 m A 



x 



T 
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2. Let us suppose that interactions CP are invariant, i.e. 



Wx^f,f' — W x ^,ff' 



In this situation from (j40|) . (|42[) follows again (f53|) . Thus and in situation of 
CP-invariance of interactions baryogenesis doesn't execute. 



3. Integrating kinetic equations for X-bosons (|42|) over momentum space and 
combining this result with equations (1401) taking into account initial conditions 
(|43[) we will receive two strict implications: 



AC 



+ ^An 2/3 + ^An a - ^An K - An e J = (54) 
law of conservation of electric charge e a n a = 0) 



E 

A,C 



2 

3 V An - 4 /3 + An -i/3 + An 2 / 3 ) 
-- (An Q + An K ) - An e - An v 



= (55) 



- law of conservation of difference of baryon and lepton charges where incorpo- 
rated notation: 

An a = n a - n a . (56) 
Summation in (|54|) . (|55p is carrying out by all charms and colors of particles. 

4- Supposing now that all X-bosons decayed i.e., proceeding to limit t — > oo and 
combining correspondences (|54|). (|55p with accounting of conditions of symmetry 
(l29l) and notations (l28l). we will receive: 



u(oo) = 2g?(oo); u(oo) = 2e(oo) =>■ d(oo) = e(oo). (57) 

From (|57f follows that on completion of decay of X-bosons for each excess 
charged lepton it will corresponds one excess katoquark and two excess anoquark, 
what subject to colors will give one excess achromic baryon to one lepton. Hence 
final baryon charge of Universe will be equal to its lepton charge: 

An b (oo) = An e (oo). (58) 

Thus for finding baryon asymmetry of Universe it is enough to define only one 
from three values: 

w(oo), e(oo), d(oo). 



5. In simplest model of baryogenesis when CP-invariance violets only in de- 
cay channels of one type of X-bosons, X_4/ 3 , from (I54|) . (|55|) strictly follows 
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correspondence in all time of process: 

e(t)=d(t), (59) 

and remain only two independent values u(t) d(t). 

Concrete numerical model of baryogenesis will be published in the next ar- 
ticle. 
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